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Abstract 



We determine mass and transverse momentum spectra of dileptons produced 
in Pb+Au (158 GeV/u) collisions within a pion annihilation model. A fit 
to the data requires simultaneous mass reduction and broadening of the in- 
medium rho propagator. The introduction of a finite pion width, as required 
within self-consistent approaches to the interacting pion gas, further improves 
the agreement with the data. 
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I. INTRODUCTION 



Since a couple of years the experimental study of dilepton production in relativistic 
heavy ion collisions is becoming more and more detailed. Besides the usual mass spectrum, 
the latest results of the CERES collaboration comprise also information on the transverse 
momenta of the produced dilcptons in Pb+Au (158 GeV/u) collisions [1,2]. 

Since several hundreds of pions are produced in such a collision, it is clear that pion 
annihilation during the hot and dense initial phase plays a dominant role for low-mass 
dilepton production, apart from the known contributions due to single meson Dalitz decays 
after freeze-out. Prom the previous analyses it is indeed well known that a major part of the 
observed dilepton "excess" can be attributed to pion-pion annihilation in the dense medium 
that is created during the collision. However, the experimental electron pair mass spectrum 
can only be explained by assuming a modification of the rho meson mass and/or width in 
the medium [3-9]. 

Apart from the knowledge of in-medium meson properties, the theoretical models require 
assumptions about the space-time development of the collision. Usually rather elaborate 
transport simulations or hydrodynamical models are used [3-5]. A much cruder, but very 
transparent way is to employ an idealized standard scenario, allowing to test easily some 
qualitative assumptions about the dilepton production mechanism. In this article we con- 
front the latest experimental results to theoretical predictions of such a "standard model" 
of the space-time evolution of an ultrarelativistic heavy ion collision, namely the Bjorken 
picture of a purely longitudinal expansion [8-10], involving an initial quark-gluon plasma 
(QGP) phase, an intermediate mixed phase, and a final hadronic (pion gas) phase. 

This model involves essentially four parameters, namely the temperatures Tj, T^, and Tf, 
specifying the initial, critical, and freeze-out temperatures; as well as the ratio r of degrees 
of freedom in the quark and hadronic phase, respectively. As in Ref. [9], wc use r = 37/3 
(two flavor quark-gluon plasma versus pion gas) and the temperatures Tj = 250 MeV, 
Tc = 160 MeV, whereas we flx Tf = 140 MeV, between the extreme cases of 120 MeV and 
150 MeV that were discussed in Ref. [9]. The precise values of these parameters are perhaps 
not so important as other fundamental assumptions of the model, such as purely longitudinal 
expansion and the idealization of pure quark and pion phases. We feel, however, that all 
these factors mainly influence the overall normalization of the resulting dilepton spectrum 
and not so much its spectral shape, which is of major interest in the following theoretical 
analysis. 



We begin with a short review of the relevant expressions. The total dilepton production 
rate due to pion annihilation in lowest order of the electromagnetic interaction can be written 

as 



II. FORMALISM 
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and the thermal pion distributions 

" exp [p ■ u{x)/T{x)] - 1 ' 

where T(a;) is the local temperature in the fluid element moving with four velocity u{x). 
The matrix element 

IM^I"^ ^ inAP3^P4) -^h'"'{pi,P2) , g =Pl +P2 =P3 +P4 (4) 

can be specified in terms of the hadronic (initial state) and leptonic (final state) tensors, 

h^''{pi,P2) = - \F^{q' f {Pi - P2nPi - P2r , (5) 
^'^^(P3,P4) = 4[p^P4"+p^p^-(p3-P4 + m,V1 , (6) 

where mi is the lepton mass and F^^ is the experimentally known pion form factor. At this 
stage the vacuum matrix element is used, without any modification due to medium effects. 
Defining 

if^^(g,x) = ^^(27r)V(pi +P2 - q) h^''{pi,P2) r{pi,x)r{P2,x) , (7) 

V(g) = lli'^'^WP^+P^ - V(P3,P4) (8) 

- {q\. q,q.) ^ , ^ ^ e{q^ - Amf) (l + ^) ^ 1 - '-f , (9) 

the differential rate for the production of dileptons with four-momentum q at the space-time 
point X can be written as 

(10) 

with 
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«[M^ - (m. + m.)T (l - (=^i±^) (l - <^^^1^) , (12) 



where — and k is the center-of-mass {q — 0) momentum of the two colliding pions 
with mass-squared mf = pj. 

The expected number of lepton pairs per collision is then: 
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X JJ^{27ry5\p, +P2-q)J d'xf-'(p,,x)r(p2,x) . (13) 

Therefore the problem reduces to the evaluation of the integral of the distribution functions 
over space time. Within the Boltzmann approximation {f^{p,x) ^ exp[—p-u{x)/T{x)], 
valid for M ^ T) these integrations can even be performed analytically [8,11], leading to 

A (271)4 67r I '^^^^l Stt 7'''''^ ^ ' 
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and 



= ^(g2 - 4m^) (^1 - ^ j , (mH^HO (16) 



/ 
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where R is the effective radius of the reaction zone, Tc is the proper time corresponding to the 
beginning of the phase transition, Mt = \/ + is the transverse mass of the dilepton, 
z = Mt/T, K{z) = (z^ + %)z^K^{z)^ and Kq and are Bessel functions. The overall 
normalization is fixed by the condition for an isentropic expansion, tiR^tT^ = KdNj^/dy, 
with K ~ 0.22 for an SUf{2) plasma. One arrives finally at the dilepton spectrum due to 
pion annihilation: 
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X (^^^^^e^Xo(^c) + [K{z^) - K{zf)]^ . (18) 



Keeping, however, the Bose distributions in Eq. (13), or in any case for a temperature 
dependent form factor, the space-time integrations have to be carried out numerically, as 
described in detail in Ref. [9] . In this case one can make use of the analytical result for the 
phase-space integral appearing in Eq. (11): 



JJ^{27r)'5\p, +P2 - q)r(pi,x)r(p2,x) = 

cosh (uj^/T) — cosh (S/T) 



Snqexp (E/T) 
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cosh {uj-/T) — cosh {5/T) 



(19) 



with 



E 2kq 
mj-ml E 

^-^^2^ (21) 

where E = q ■ u{x) is the dilepton energy, and q = \JE'^ — is the three-momentum 
component of q in the production fluid element moving with four-velocity u. 

The previous expressions concern dileptons stemming from the pion phase. One has to 
add the contribution due to quark- antiquark annihilation that we have not written down 
explicitly [8,9,11], because it is not concerned by the considerations in this paper. In order to 
compare with the experimental mass and transverse momentum spectra, the sum of the two 
parts has then to be multiplied with the relevant experimental acceptance function A{M, qx) 
and folded with the experimental mass resolution, see Ref. [9]. Finally the experimentally 
determined contribution due to decays of mesons produced during the primary nucleon- 
nucleon coUisions ("cocktail") has to be added. We discuss now the results obtained after 
this procedure. 



III. RESULTS 

A first glance at the experimental data (markers in Fig. 1) in comparison to the expec- 
tation from primary meson decays (dotted line) shows that practically all excess dileptons 
that are produced in the reaction zone are required to fall within the mass region [200,600] 
MeV. In particular close to the bare rho mass (?=a 800 MeV) practically no additional contri- 
bution is required. This fact indicates already a substantial modification of the in-medium 
pion form factor, which is indeed confirmed by our theoretical result using the vacuum pion 
form factor (solid curve), which largely overshoots the experimental results in the rho mass 
region. An evaluation using the pure Born cross section, i.e., setting tentatively F^^ — 1 
(dot-dashed curve) yields instead a fair agreement in the rho mass region with, however, 
missing strength at lower mass. It is therefore clear that necessarily a substantial reduction 
of the rho meson mass in the medium is required in order to reproduce the experimental 
data. 

Since in our model all excess dileptons stem from a pion gas phase with a small vari- 
ation of temperature between Tc = 160 and Tf = 140 MeV, we neglect even the detailed 
dependence on temperature and rather attempt to construct an effective in-medium pion 
form factor of the standard functional form 

\FAM)f = J—, (22) 

that allows to reproduce the experimental data. The constants nip, Fp, and C are now 
considered as free parameters and allowed to deviate from their vacuum values m° = 0.76 
GeV, = 0.135 GeV, and C° = 1.3. Indeed we find that a choice of = OAm% Tp = 2T% 
and C = 7C° allows a very reasonable fit of the data, shown in Fig. 2 (solid line) . Thus a 
simultaneous modification of centroid (reduction by a factor 2 ... 3) and width (increase by 
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about a factor 2) of the pion form factor is required; we found that changing only one of 
them does not yield satisfactory results. The comparison of the resulting effective in-medium 
and the vacuum pion form factor is shown in Fig. 3. We remark at this point that in the 
present model the contribution of dileptons from the quark-gluon plasma (produced during 
the initial phase and the mixed phase) is neghgible [7,8,12,13], as shown by the dashed hne 
in Figs. 1 and 2. This is mainly due to the short duration of this phase compared to the 
pion phase. 

After discussing the mass spectrum of electron pairs, we come now to the distribution 
of transverse momenta. Fig. 4 shows the experimental values (in the mass interval of inter- 
est 250 MeV < M < 680 MeV) in comparison with our theoretical prediction (solid line). 
The good agreement is indeed remarkable, given the fact that in our model the transverse 
momentum of the dileptons is entirely due to the thermal motion of the colliding pions in 
the longitudinally expanding medium. This is the reason for the larger slope of the trans- 
verse momentum distribution of dileptons stemming from pion annhihilation (dashed line) 
compared to the Dahtz dileptons (dotted line) . The sum of both results in fair agreement 
with the data. It seems thus that the transverse momenta of the dileptons reveal no hint of 
a transverse expansion of the reaction zone. 

The in-medium pion annihilation contribution to the dilepton mass spectrum is naturally 
cut off at the two-pion threshold, and leaves therefore a small dip in the theoretical curves 
around M 300 MeV (sohd lines in Figs. 1 and 2). A natural remedy would be a decrease 
of the pion mass in the medium, for which, however, there is so far no solid theoretical 
foundation. In this article we rather study the consequences of introducing empirically a 
finite damping width of the pion, due to resonant pion rescattering in the medium. The 
occurence of a finite pion damping width F^r ~ 50 MeV is well known from the study of 
the ttA^A system [14] and is also expected to occur in a hot dense pion gas, due to the 
strong coupling to the rho and ai resonances. Within a simple field theoretical model for 
this system, motivated by a quartic self-coupling of the pions, van Hees and Knoll [15] have 
shown how the self-consistent equations of motion can be derived using Baym's $ functional 
approach [16]. As a result of these models, a temperature and momentum dependent pion 
width has been obtained [17], which in the relevant domain of T 150 MeV amounts to 
Ftt ~ 50 . . . 100 MeV. We therefore assume a Gaussian pion mass distribution of width Ftt, 



9n{m) = exp 
v27ri „ 



(m — m,r)^ 



2r2 



(23) 



and naively generalize Eq. (14) to the reaction of two pions with different masses 7r(mi) -|- 
7r(m2) — >■ 7(?)- This is achieved by replacing the function i?(g^), Eq. (16), by 



/•oo /"OO 

Hr„{q^) = / dmig^{mi) / dm2 gn{m2) eiq"^ - {mi + 1712)^] 
Jo Jo 



X 



^ _ {mi+m2y \ / (mi - 1712^ 



3/2 



(24) 



which is plotted in Fig. 5 for two choices of the pion width F^r = 50, 100 MeV. The resulting 
dilepton spectrum for a value F^ = 100 MeV (and form factor parameters nip — rrip/S, 
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Tp = 2r°, and C — 40C'') is also shown in Fig. 2 (solid line). Indeed the introduction of a 
finite pion width can perfectly describe the spectrum also in the threshold region M ^ 2m^. 

Of course, there are other processes that could potentially contribute dileptons in this 
mass region around 300 MeV. In particular, we neglected completely the possible presence 
of baryons in the reaction zone and of baryon induced reactions like pN — > B, which might 
contribute in the low-mass region. Also possible mesonic rescattering of the rho resonance, 
pTT — > 7re+e_, has been disregarded [18]. However, since the mass spectrum below 200 MeV 
is perfectly reproduced by ttq, t], and r]' Dalitz decays [19], there is not much freedom for 
such additional "exotic" contributions, since these would dominantly contribute in this low- 
mass domain, which is on the contrary not the case for the pion width mechanism discussed 
above. Indeed also theoretical estimates of these processes predict too small rates [12,20]. 

IV. CONCLUSIONS 

In summary, the latest Pb+Au CERES data hint to a simultaneous reduction of the rho 
meson mass and an increase of its width in the dilepton production medium. Even if our 
model (Bjorken scenario) of the space-time evolution of the collision is rather idealized, we 
think that these qualitative conclusions based on the shape of the dilepton mass distribution 
can be drawn. The further ad-hoc introduction of a finite mass distribution of the interacting 
pions served to simulate the effect of a finite pion damping width, which is required for a 
consistent treatment of rho mesons as decaying resonances in an interacting pion gas. It 
allowed us to reproduce the experimental results around the two-pion threshold even better. 

Surprisingly, the transverse momenta of the dileptons were also rather well described 
by our model, and thus appear to reveal no hint on a transverse expansion of the reaction 
zone. Finally, dilepton production from an eventual QGP phase is negligible due to the short 
duration of such a phase, even if its existence is an assumption of the chosen theoretical model 
and influences indirectly, through a consistent time evolution, also the hadronic rates. 

The present investigation can be improved by applying a more realistic hydrodynamical 
evolution which goes beyond the simple Bjorken scenario. The pion broadening should be 
considered within a definite microscopic model, so that its temperature and momentum 
dependence can be included into the calculation. Such a model should allow to calculate the 
pion broadening consistently together with the rho meson spectral properties. Since already 
the simple calculation of the rho meson spectral function with on-shell pions in a hot pion 
gas yields a downwards shift of the rho meson mass by about 30 MeV [21], it can be expected 
that off-shell pions with about 100 MeV width could explain an even larger rho meson mass 
shift, which is inevitable in order to explain the new low- mass dilepton data of the CERES 
collaboration. The objective of this article was, however, the empirical determination of 
the characteristic magnitudes of the in-medium pion and rho meson parameters that should 
result from a theoretical model in order to explain the present experimental data. 

The findings of our analysis may be carried over to a prediction of dilepton spectra for 
the conditions at RHIC, where our pion annihilation model together with the QCD phase 
transition and the suggested in-medium modifications of pion and rho meson properties will 
be appropriate. 
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FIGURES 
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FIG. 1. The dilepton mass spectrum evaluated with a free-space (solid line) or a Born 

(dot-dashed line) pion form factor. The theoretical predictions involve acceptance and resolu- 
tion corrections, and the primary meson decay contribution (dotted line) is added. The dashed 
line shows separately the contribution of qq annihilation in the quark and mixed phase. The data 
points and the meson decay contribution are taken from Ref. [2]. 
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FIG. 2. As Fig. 1, but with the in- medium pion form factor (solid line), and additionally 
assuming a finite pion width (dot-dashed line). 
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FIG. 3. Free-space and in-medium pion form factor. 
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FIG. 4. The dilepton transverse momentum spectrum. 
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FIG. 5. Kinematic factor, Eq. (24). 



11 



